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Abstract 

We use the density matrix renormalization group (DMRG) method to study the 
ground and low-lying excited states of three kinds of uniform and dimerized alternat- 
ing spin chains. The DMRG procedure is also employed to obtain low-temperature 
thermodynamic properties of these systems. We consider a 2N site system with 
spins si and S2 alternating from site to site and interacting via a Heisenberg an- 
tiferromagnetic exchange. The three systems studied correspond to (si,S2) being 
equal to (1, (§, 5) an d (§ 5 1); all of them have very similar properties. The 
ground state is found to be ferrimagnetic with total spin sq = N(si — ^2). We 
find that there is a gapless excitation to a state with spin sq — 1, and a gapped 
excitation to a state with spin sq + 1- Surprisingly, the correlation length in the 
ground state is found to be very small for this gapless system. The DMRG analysis 
shows that the chain is susceptible to a conditional spin-Peierls instability. Further- 
more, our studies of the magnetization, magnetic susceptibility x an d specific heat 
show strong magnetic-field dependences. The product xT shows a minimum as a 
function of temperature T at low magnetic fields; the minimum vanishes at high 
magnetic fields. This low-field behavior is in agreement with earlier experimental 
observations. The specific heat shows a maximum as a function of temperature, and 
the height of the maximum increases sharply at high magnetic fields. Although all 
the three systems show qualitatively similar behavior, there are some notable quan- 
titative differences between the systems in which the site spin difference, |si — S2I, 
is large and small respectively. 

PACS number: 75.50. Gg 
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1 Introduction 



There have been a number of experimental efforts over the last several years to synthesize 
molecular systems showing spontaneous magnetization 0, H Hi- These studies have led 
to the recent realization of systems such as NPNN (para-nitrophenyl nitronyl nitrox- 
ide) and 6*60 — TDAE (Tetrakis Dimethyl Amino Ethylene) which show ferromagnetic to 
paramagnetic transition at low-temperatures. There also exist theoretical models |IJ [IJ 
that predict ferrimagnetism in organic polymer systems, but a truly extended organic 
ferrimagnet has not been synthesized although an oligomer with a ground state spin of 
S = 9 is now known. On the inorganic front, pursuit of molecular magnetism has been 
vigorous, and there has been success in the synthesis of molecular systems showing sponta- 
neous magnetization at low-temperatures || [7j. These are quasi-one-dimensional bimetal- 
lic molecular magnets in which each unit cell contains two spins with different spin values 
0. These systems contain two transition metal ions per unit cell, with the general for- 
mula ACu(pbaOH)(H20)3.2H20, where pbaOH is 2-hydroxo-l,3-propylenebis(oxamato) 
and A=Mn, Fe, Co, Ni, and they belong to the alternating or mixed spin chain family 
||. These alternating spin compounds have been seen to exhibit ferrimagnetic behavior. 
It has been possible to vary the spin at each site from low values where quantum effects 
dominate to large values which are almost classical. 

There have been a number of theoretical investigations for quantum ferrimagnetic 
systems in recent years |9|, [10|. However, these models consider complicated multispin 



interactions, while very little is known about the simple Heisenberg model with purely 
quadratic interactions. There have been some analytic studies of alternating spin-l/spin-i 
system [[T| and a detailed spin-wave analysis followed by a DMRG study [0 correspond- 
ing to the Cu — Ni bimetallic chain with the simple Heisenberg model. These studies 
confirm that these ferrimagnetic systems can be accurately described by a pure Heisenberg 
spin model. 

Theoretical studies of alternating spin systems have so far been concerned only about 
the spin- 1 /spin- i bimetallic chains; further, the thermodynamic properties have not been 
explored in detail. The thermodynamic behavior of these alternating spin compounds is 



very interesting || [13]. In very low magnetic fields, these systems show one-dimensional 
ferrimagnetic behavior. The yT vs. T (where \ is the magnetic susceptibility and T 
the temperature) plots show a rounded minimum. As the temperature is increased, yT 
decreases sharply, goes through a minimum before increasing gradually. The temperature 
at which this minimum occurs differs from system to system and depends on the site spins 
of the chain. The variation of the field induced magnetization with temperature is also 
interesting as the ground state is a magnetic state. At low magnetic fields, the magneti- 
zation decreases with increasing temperature at low temperatures. However, at moderate 
magnetic fields, with increase in temperature, the magnetization slowly increases, shows 
a broad peak and then decreases. Such a behavior has been studied theoretically by us 
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for the spin-l/spin-i system flT2 



These interesting observations have motivated us to study of ferrimagentic systems 
with arbitrary spins Si and S2 alternating from site to site. It would be quite interesting 
to know the thermodynamic properties of these systems with varying s\ and $2- In the 
previous paper, we predicted some interesting features of the thermodynamic properties 
of the spin- 1 /spin- i system at high magnetic fields. We focus on the high field behavior 
of the general ferrimagnetic chains to find whether the observed properties are generic to 
these systems or are dependent on specific s\ and S2 values. 

In this paper, we study the low-lying excited states and low-temperature properties 
of the spin- 1 /spin- 1 and spin- 1 /spin- 1 chains and rings and compare them with with 
those of the corresponding spin-l/spin-| systems. We have employed the density matrix 
renormalization group (DMRG) method which has proved to be the best numerical tool 
for low- dimensional spin systems in recent years jnj. The ground state energy per site, 
the spin excitation gap and the two-spin correlation functions obtained from this method 
have been found to be accurate to several decimal places [fl5| , [16], |I7]] when compared with 
Bethe-ansatz results (where possible) and exact diagonalization results of small systems. 
In the DMRG method, spin parity symmetry can be used to characterize the spin states 
along with the S^ ot as the good quantum numbers. The DMRG calculations have been 
carried out on chains and rings with alternate spin-si/spin-S2 sites, with s\ fixed at | and 
s 2 being 1 or ~. Studies of the ground state and low-lying excited states are reported in 
detail and compared with those of Si = 1 and S2 = \ system. Furthermore, by resorting 
to a full diagonalization of the DMRG Hamiltonian matrix in different S z sectors, we 
have also obtained the low-temperature thermodynamic properties of these systems. The 
theromo dynamic properties we discuss will include low and high field magnetization, 
magnetic susceptibilty and specific heat, all at low temperatures. 

The paper is organized as follows. In section 2, we present properties of the ground 
and low-lying excited states and compare them with the results of spin wave theory. In 
section 3, we discuss the low-temperature thermodynamic properties of the systems. We 
summarize our results in the last section. 



2 Ground State and Excitation Spectrum 

We start our discussion with the Hamiltonian for a chain with spins Si and S2 on alter- 
nating sites (with s\ > S2, without loss of generality). 

H = J ^2 [(1 + 8) S l!n ■ S 2 , n + (1-5) S 2 , n ■ S hn+1 } , (1) 

n 

where the total number of sites is 2N and the sum is over the total number of unit cells 
N. Si tn corresponds to the spin operator for the site spin Si in the n-th unit cell. The 
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exchange integral J is taken to be positive for all our calculations; 5 is the dimerization 
parameter and it lies in the range [0, 1]. 

2.1 Summary of results from spin- wave theory 

Before describing our numerical results, we briefly summarize the results of a spin wave 
analysis for the purposes of comparison jn|. We will first state the results with 5 = 0. 
According to spin wave theory, the ground state has total spin sq = N(si — s 2 ). Let us 
define a function 

uj(k) = J ^/(si -s 2 ) 2 + 4s lS2 sin 2 (fc/2) , (2) 
where k denotes the wave number. Then the ground state energy per site is given by 

En 1 r dk , , , . , . . 

e = = - Js lS2 + - — [ - J(s 1 + s 2 + u{k) ] . (3) 

2iv Z Jo ix 

The lowest branch of excitations is to states with spin s = sq — 1, with the dispersion 

uJi(k) = J(-s 1 + s 2 ) + u(k) ; (4) 

the gap vanishes at k = 0. There is a gapped branch of excitations to states with spin 
s = sg + 1, with the dispersion 

uj 2 {k) = J (si - s 2 ) + ; (5) 

the minimum gap occurs at k = and is given by A = 2J(s\ — s 2 ). In the ground state 
with S z = Sq, the sublattice magnetizations are given by the expectation values 

1 1 r* dk J(si + s 2 ) 

s i + o ) ~ o / W ' 

2 z Jo i co> 

Sl - S2 - < S z hn > . (6) 

The various two-spin correlation functions decay exponentially with distance; the inverse 
correlation length is given by £ -1 = ln(si/s 2 ). The results with dimerization (5 > 0) are 
very similar. In fact, within spin wave theory, the minimum gap A to states with spin 
s = sg + 1 is independent of S. 

2.2 Results from DMRG studies 

We have studied the system defined by equation ([J) both with and without dimerization, 
5^0 and 5 = respectively. We study alternating spin-|/spin-l (hereafter designated 
as (|,1)) and spin-|/spin-| (to be called (|, |)) chains with open boundary condition 
for the Hamiltonian (fj) by employing the DMRG method. We compute the ground 
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state properties for both the systems by studying chains of upto 80 to 100 sites. The 
number of dominant density matrix eigenstates, m, that we have retained at each DMRG 
iteration also varies between 80 to 100 for both the systems. With the increase of the Fock 
space dimensionality of the site spins, we increase m to obtain more accurate results. The 
DMRG procedure follows the usual steps for chains discussed in earlier papers ]I4|, [L6|, [17|] , 
except that the alternating chains studied here are not symmetric between the left and 
right halves; hence the density matrices for these two halves have to be constructed at 
every iteration of the calculations. We have verified the convergence of our results by 
varying the values of m and the system size. The ground states of both the system lie in 
the S z = N(s\ — S2) sector, as verified from extensive checks carried out by obtaining the 
low-energy eigenstates in different S z sectors of a 20-site chain. A state corresponding to 
the lowest energy in S z = N(si — S2) is found in all subspaces with \S Z \ < N(si — S2), and 
is absent in subspaces with \S Z \ > N(si — S2). This shows that the spin in the ground 
state is sg = N(si — S2). 

In a previous paper, we studied the alternating spin chain made up of spin-1 and 
spin- 1 (hereafter refered as (1, 5)) a t alternate sites flT2| . We compare here the results for 
all the three systems, nemely, (|, 1), (|, 1/2) and (1, ~) spin systems. In figure 1, we show 
the extrapolation of the energy per site as a function of inverse system size for all three 
systems. The ground state energy per site e extrapolates to — 1.93096J for the (|, 1) 
system, to -0.98362J for (§ A), and to -0.72704J for (1, |). The spin wave analysis 
gives the values -1.9147 for (§, 1), -0.979J for (§, §) and -0.718J for (1, \) which are 
all higher than the DMRG values. It is interesting to note that, in the alternating spin 
cases, the energy per site lies in between the values for the pure spin-si uniform chain 
and the pure spin-S2 uniform chain. 

In figure 2, we show the expectation value of site spin operator Sf n ( spin density ) 
at all the sites for the (|, 1), (§, |) and (1, §) chains. The spin densities are uniform on 
each of the sublattices in the chain for all the three systems. For the (|, 1) chain, the spin 
density at a spin-| is 1.14427 (the classical value is |), while, at a spin-1 site it is —0.64427 
(classical value 1). For the (|, |) case, the spin density at a spin-| site is 1.35742 and at 
a spin-A site it is —0.35742. For the (1, A) case, the value at a spin-1 site is 0.79248 and 
at a spin-| site it is —0.29248. These can be compared with the spin wave values of 1.040 
and —0.540; 1.314 and —0.314; and 0.695 and —0.195 for the spin-si and spin-s 2 sites of 
the (|,1); (|, |); and (1, ~) systems respectively. We note that the spin wave analysis 
overestimates the quantum fluctuations in case of systems with small site spin values. We 
also notice that there is a greater quantum fluctuation when the difference in site spin 
|«i — s 2 | is larger. This is also seen in spin- wave theory. The spin density distribution in an 
alternating (si, S2) chain behaves more like that in a ferromagnetic chain rather than like 
an antiferromagnet, with the net spin of each unit cell perfectly aligned (but with small 
quantum fluctuations on the individual sublattices). In a ferromagnetic ground state, the 
spin density at each site has the classical value appropriate to the site spin, whereas for an 
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antiferromagnet, this averages out to zero at each site as the ground state is nonmagnetic. 
From this viewpoint, the ferrimagnet is similar to a ferro magnet and is quite unlike an 
antiferromagnet. The spin wave analysis also yields the same physical picture. 

Because of the alternation of spin-si and spin-s 2 sites along the chain, one has to 
distinguish between three different types of pair correlations, namely, < S^ S^ n >, < 
S20S2 n > and < Sl fi S^ n >• We calculate all the three correlation functions with the mean 
values subtracted out, since the mean values are nonzero in all these three systems unlike 
in pure antiferromagnetic spin chains. In the DMRG procedure, we have computed these 
correlation functions from the sites inserted at the last iteration, to minimize numerical 
errors. In figure 3, we plot the two-spin correlation functions in the ground state as a 
function of the distance between the spins for an open chain of 100 sites for all three 
cases. All three correlation functions decay rapidly with distance for each of the three 
systems. From the figure it is clear that, except for the < Sf 0^2 n > correlation, all other 
correlations are almost zero even for the shortest possible distances. The < S^S^n > 
correlation has an appreciable value [—0.2 for (|,1), —0.07 for (§, |) and —0.094 for 
(1, i)] only for the nearest neighbors. This rapid decay of the correlation functions makes 
it difficult to find the exact correlation length £ for a lattice model, although it is clear 
that £ is very small (less than one unit cell) for the (|, |) and (1, |) cases, and a little 
greater (1 < £ < 2) for the (|, 1) system. Spin wave theory gives £ = 2.47 for (|,1), 
£ = 0.91 for (§, I), and £ = 1.44 for (1, \) cases. 

The lowest spin excitation of all the three chains is to a state with s — sq — 1- To 
study this state, we target the 2 nd state in the S z = s G — l sector of the chain. To confirm 
that this state is a s = sq — 1 state, we have computed the 2-nd state in S z — sector 
and find that it also has the same energy. However, the corresponding state is absent 
in S z sectors with \S Z \ > s G — 1. Besides, from exact diagonalization of all the states 
of all the si — s 2 alternating spin chains with 8 sites, we find that the energy orderings 
of the states is such that the lowest excitation is to a state with spin s — sq — 1- We 
have obtained the excitation gaps for all the three alternating spin chains in the limit 
of infinite chain length by extrapolating from the plot of spin gap vs. the inverse of the 
chain length (figure 4). We find that this excitation is gapless in the infinite chain limit 
for all three cases. This property is quite unusual; the spin systems studied till date with 
antiferromagnetic exchage interactions have a small correlation length if and only if there 
is a finite gap to the lowest excited state of the system, from the ground state. Thus, this 
gapless excitation coexisting with very short correlation length in the ground state can 
indeed be taken as a signature of systems with spontaneous magnetization. 

To characterize the lowest spin excitations completely, we also have computed the 
energy of the s = sq + 1 state by targetting the lowest state in the S z — sq + 1 sector. 
In figure 5, we plot the excitation gaps to the s = sq + 1 state from the ground state for 
all three systems as a function of the inverse of the chain length. The gap saturates to 
a finite value of (1.0221 ± 0.0001) J for the (§, 1) case, (1.8558 ± 0.0001) J for (§, §), and 
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(1.2795 ± 0.0001) J for (1, |). It appears that the gap is also higher when the difference in 
site spins \s\ — S2I, is larger. The site spin densities expectation values computed in this 
state for all three cases are found to be uniform (i.e., independent of the site) on each of 
the sublattices. This leads us to believe that this excitation cannot be characterized as 
the states of a magnon confined in a box, as has been observed for a spin-1 chain in the 
Haldane phase [|l^ . 



Earlier studies on pure spin-1 and pure spin-| chains |L8] have revealed that with 
the alternation 5 in the exchange parameter, the half-odd-integer spin chain will have 
an unconditional spin-Peierls transition while for integer spin chains, the transition is 
conditional. This conclusion has been drawn from the fact that, with the inclusion of 5, 
the magnetic energy gain AE can be defined as 

AE(2N, 5) = ^ [ E(2N, 5) - E(2N, 0) ] , (7) 

where E(2N, 5) is the ground state energy of the 2iV-site system with an alternation S 
in the exchange integral, and E(2N, 0) is the ground state energy of the uniform chain 
of 2N sites. For the pure spin chain, if we assume that AE varies as S u for small 5, we 
find that v < 2 for the half-odd-integer spin chains and v = 2 for the integer spin chains 
Thus, for a half-odd-integer spin chain, the stabilization energy always overcomes 



the elastic energy, whereas for the integer spin case, it depends on the lattice stiffness. 

We have used DMRG calculations extend these studies to the alternating spin systems. 
We obtain AE(2N, S), for small values of 5 for all the three alternating spin chains. To 
determine the exact functional form of the magnetic energy gain, we varied the chain 
length from 50 sites to 100 sites and also m (the number of states retained in each DMRG 
iteration) from 80 to 100 to check the convergence of AE with chain length. We note that 
as we approach the classical spin limit, the convergence is reached faster. The dependence 
of AE(2N, S) on 1/2N is linear for all three cases for the S values we have studied. Figure 
6 gives a sample variation of AE(2N, S) with 1/2N for all the systems we have studied. 
This allows us to extrapolate AE(2N, 5) to the infinite chain limit reliably for all the 
cases. In figure 7, we show the plot of AE(2N, 5) vs. S for finite 2N values and also the 
extrapolated infinite chain values for all three systems. We see that there is a gain in 
magnetic energy upon dimerization even in the infinite chain limit for all the systems. To 
obtain the exponent u, we plot In AE(2N, 5) vs. In 5 for the infinite chain (figure 8). From 
this figure, we find that in the alternating spin case, for the infinite chain AE « $2.oo±o.oi 
for all three cases. Thus the spin-Peierls transition appears to be close to being conditional 
in these systems. The magnetic energy gain per site for finite chains is larger than that 
of the infinite chain for any value of S (figure 7). It is possible that the distortion in a 
finite chain is unconditional while that of the infinite chain is conditional for ferrimagnetic 
systems. 

We have also studied the spin excitations in the dimerized alternating (si,S2) chains, 
defined in equation ([!]). We calculate the lowest spin excitation to the s = s^ — 1 state from 
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the ground state. We find that the s = sq — 1 state is gapless from the ground state for all 
values of 5. This result agrees with the spin wave analysis of the general (si, S2) chain. The 
systems remain gapless even while dimerized unlike the pure antiferromagnetic dimerized 
spin chains. There is a smooth increase of the spin excitation gap from the ground state 
to the s = sq + 1 state with increasing 5 for all three systems studied here. We have 
plotted this gap vs. 5 in figure 9. The gap shows almost a linear behavior as a function of 
S, with an exponent of 1.0±0.01 for all three systems |I9"|| . This seems to be an interesting 
feature of all ferrimagnets. The spin wave analysis however shows that this excitation gap 
is independent of 5 for the general (si,S2) chain. The similar behaviors of these three 
alternating spin systems suggest that a ferrimagnet can be considered as a ferromagnet 
with small quantum fluctuations. 



3 Low- Temperature Properties 

In this section, we present results of our DMRG calculations of the thermodynamic prop- 
erties of the (|,1), (§, |) and (1, ~) spin systems. The size of the system varies from 8 
to 20 sites. We impose periodic boundary conditions to minimize finite size effects with 
Si,n+i = so that the number of sites equals the number of bonds. We set up the 
Hamiltonian matrices in the DMRG basis for all allowed S z sectors for a ring of 2N sites. 
We can diagonalize these matrices completely to obtain all the eigenvalues in each of the 
S z sectors. As the number of DMRG basis states increases rapidly with increasing m, we 
retain a smaller number of dominant density matrix eigenvectors in the DMRG proce- 
dure, i.e., 50 < m < 65, depending on the S z sector as well as the size of the system. We 
have checked the dependence of properties (with m in the range 50 < m < 65) for the 
system sizes we have studied (8 < 2N < 20), and have confirmed that the properties do 
not vary significantly for the temperatures at which they are computed; this is true for 
all the three systems. The above extension of the DMRG procedure is found to be ac- 
curate by comparing with exact diagonalization results for small systems. It may appear 
surprising that the DMRG technique which essentially targets a single state, usually the 
lowest energy state in a chosen sector, should provide accurate thermodynamic properties 
since these properties are governed by energy level spacings and not the absolute energy 
of the ground state. However, there are two reasons why the DMRG procedure yields rea- 
sonable thermodynamic properties. Firstly, the projection of the low-lying excited state 
eigenfunctions on the DMRG space in which the ground state is obtained is substantial 
and hence these excited states are well described in the chosen DMRG space. The second 
reason is that the low-lying excitations of the full system are often the lowest energy states 
in different sectors in the DMRG procedure; hence their energies are quite accurate even 
on an absolute scale. 
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The canonical partition function Z for the 2N site ring can be written as 



Z = J2 e-ft E *- B Mi) , (8) 

3 

where the sum is over all the DMRG energy levels of the 2N site system in all the S z 
sectors. Ej and (M)j denote the energy and the ^-component of the total spin of the state 
j, B is the strength of the magnetic field in units of J/gns {g is the gyromagnetic ratio 
and fx b is the Bohr magneton) along the z direction, and (3 = J/ksT with ks and T being 
the Boltzmann constant and temperature respectively. The field induced magnetization 
< M > is defined as 

< M > = U [M)j \ . (9) 

Z 

The magnetic susceptibility x is related to the fluctuation in magnetization 

X = PI <M 2 >-<M> 2 ] . (10) 

Similarly the specific heat C is related to the fluctuation in the energy and can be written 

as 

C = ^[<E 2 >-<E> 2 }. (11) 

The dimensionalities of the DMRG Hamiltonian matrices that we completely diag- 
onalize vary from 3000 to 4000, depending upon the DMRG parameter m and the S z 
value of the targetted sector, for rings of sizes greater than 12. These matrices are not 
very sparse, owing to the cyclic boundary condition imposed on the system. The DMRG 
properties compare very well with exact results for small system sizes amenable to exact 
diagonalization studies. In the discussion to follow, we present results on the 20-site ring 
although all calculations have been carried out for system sizes from 8 to 20 sites. This is 
because the qualitative behavior of the properties we have studied are similar for all the 
ring sizes in this range for all three systems. 

We present the dependence of magnetization on temperature for different magnetic 
field strengths in figure 10 for all three systems. At low magnetic fields, the magnetization 
shows a sharp decrease at low temperatures and shows paramagnetic behavior at high 
temperatures. As the field strength is increased, the magnetization shows a slower decrease 
with temperature, and for high field strengths the magnetization shows a broad maximum. 
This behavior can be understood from the type of spin excitations present in these systems. 
The lowest energy excitation at low magnetic fields is to a state with spin s less than sq- 
Therefore, the magnetization initially decreases at low temperatures. As the field strength 
is increased, the gap to spin states with s > s G decreases as the Zeeman coupling to these 
states is stronger than to the states with s < sq- The critical field strengths at which 
the magnetization increases with temperature varies from system to system since this 
corresponds to the lowest spin gap of the corresponding system. The behavior of the 
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system at even stronger fields turns out to be remarkable. The magnetization in the 
ground state (T = 0) shows an abrupt increase signalling that the ground state at this 
field strength has S z > sq- The temperature dependence of the magnetization shows a 
broad maximum indicating the presence of states with even higher spin values lying above 
the ground state in the presence of this strong field. In all three cases, the ground state 
at very high field strengths should be ferromagnetic. For the systems at such high fields, 
the magnetization decreases slowly with increase of temperature as no other higher spin 
states lie above the ground state. While, we have not studied such high field behaviors, 
we find that the field strength corresponding to switching the spin of the ground state sc 
to sg + 1 is higher for (|, \) system compared to (§, 1) and (1, \) systems. The switching 
field appears to depend on the value of \s\ — «2 1 ■ We see in figure 10 that in the (|, 1) and 
(1, |) cases, the ground state has switched to the higher spin state at highest magnetic 
field strength we have studied but in the (|, |) case, the ground state has not switched 
even at the field strength indicating that the excitation gap for this system is the larger 
than the other two. For the (|, \) case, the same situation should occur at very high 
magnetic fields. Thus, we predict that the highest s z is attained in the ground state at 
high magnetic field and this field strength increases with increase in site spin difference 
|si - s 2 \. 

The dependence of magnetization on the magnetic field is shown at different tempera- 
tures in figure 11 for all three systems studied here. At low temperature the magnetization 
shows a plateau. The width of the plateau depends on the system and it decreases as the 
temperature is raised. Eventually the plateau disappears at higher temperatures. The 
existence of the plateau shows that the higher spin states are not accessible at the chosen 
temperature. The plateau is widest in the (|, |) case once again reflecting the larger gap 
to the sg + 1 state in this system compared with the (| — 1) and (1, |) systems. At higher 
fields, the larger Zeeman splittings of higher spin states become accessible leading to an 
increase in the magnetization. The magetization curves at all temperatures intersect at 
some field strengths depending on the system. For the (|, 1) case, the intersection occurs 
at B = 1.0 J / g^B an d at higher field strength all the curve collapse. For (|, |) case, the 
intersection occurs at B — 1.5 J/ gfis, while for the (1, |) system these curves intersect 
twice, once at B = 1.0J/gfi B and again at B = 2.5J/gfis in the chosen field range. These 
fields are close to the field strengths at which the ground state switches from one S z value 
to a higher value for the corresponding system. Thus, they are strongly dependent upon 
the actual values of Si and S2- 

The dependence of xl '/2N on temperature for different field strengths are shown 
in figure 12 for all three systems. For zero field, the zero temperature value of xT is 
infinite in the thermodynamic limit; for finite rings it is finite and equal to the average 
of the square of the magnetization in the ground state. For the ferrimagnetic ground 
state \T jlN , as T — > 0, is given by sg{sg + l)/6iV. As the temperature increases, this 
product decreases and shows a minimum before increasing again. For the three systems 
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studied here, the minimum occurs at different temperatures depending on the system. 
For the (|, 1) alternating spin system, it is at k B T = (0.8 ±0.1) J, while for the (|, \) and 
(1, \) cases, it occurs at k B T = (1.0 ± 0.1)J and k B T = (0.5 ± 0.1) J respectively. The 
minimum occurs due to the states with S z < sc getting populated at low temperatures. 
In the infinite chain limit, these states turn out to be the gapless excitations of the 
system. The subsequent increase in the product yT is due to the higher energy-higher 
spin states being accessed with further increase in temperature. This increase is slow in 
(|, |) case, as in this system very high spin states are not accessible within the chosen 
temperature range. Experimentally, it has been found in the bimetallic chain compounds 
that the temperature at which the minimum occurs in the yT product depends upon the 
magnitude of the spins s\ and s% [|I3"|| . The Ni 11 — Cu 11 bimetallic chain shows a minimum 
in yT/2N at a temperature corresponding to 55 cm -1 (80i^); an independent estimate 
of the exchange constant in this system is 100 cmT x |[20|| . This is in very good agreement 
with the minimum theoretically found at temperature (0.5 ± 0.1) J for the (1, |) case. 
The minimum in yT/2N vanishes at B = 0.1J/g{i B which corresponds to about 10T for 
all three systems. It would be interesting to study the magnetic susceptibility of these 
systems experimentally under such high fields. The low-temperature zero-field behavior 
of our systems can be compared with the one-dimensional ferromagnet. In the latter, the 



spin wave analysis shows that the yT product increases as 1/T at low temperatures [21 

In finite but weak fields, the behavior of yT is different. The magnetic field opens up 
a gap and yT falls exponentially to zero for temperatures less than the gap in the applied 
field for all three systems. Even in this case a minimum is found at the same temperature 
as in the zero-field case for the corresponding system, for the same reason as discussed in 
the zero field case. 

In stronger magnetic fields, the behavior of yT from zero temperature upto k B T = 
Jmin (Jmin is the temperature at which the minimum in yT is observed) is qualitatively 
different. The minimum in this case vanishes for all three systems. In these field strengths, 
the states with higher S z values are accessed even below k B T = J m i n . The dependence 
of yT above k B T = J min at all field strengths is the same in all three systems. In even 
stronger magnetic fields, the initial sharp increase is suppressed. At very low temperature, 
the product yT is nearly zero and increases almost linearly with T over the temperature 
range we have studied. This can be attributed to a switch in the ground state at this 
field strength. The very high temperature behavior of yT should be independent of field 
strength and should saturate to the Curie law value corresponding to the mean of magnetic 
moments due to spin-si and spin-s 2 . 

The temperature dependence of specific heat also shows a marked dependence on the 
magnetic field at strong fields. This dependence is shown in figure 13 for various field 
strengths for all the three systems. In zero and weak magnetic fields, the specific heat 
shows a broad maximum at different temperatures which are specific to the system. Inter- 
estingly, the temperature at which the specific heat shows a maximum closely corresponds 
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to the temperature where the low-field yT has a minimum for the corresponding system. 
For a strong magnetic field (B = J), there is a dramatic increase in the peak height at 
about the same temperature corresponding to the specific system, although the qualitative 
dependence is still the same as at low magnetic fields in all three cases. This phenomena 
indicates that the higher energy high-spin states are brought to within ksT of the ground 
state at this magnetic field strength for all three cases. 

Studies of thermodynamic properties of the dimerized alternating spin chains in these 
three cases show qualitatively similar trends to that of the corresponding uniform systems 
and follows from the fact that the low-energy spectrum does not change qualitatively upon 
dimerization. 

4 Summary 

We have studied alternating spin-|/spin-l and spin-|/spin-i systems in detail and com- 
pared them with our earlier studies on spin-l/spin-| model. The ground state and low- 
lying excited states have been analyzed by employing the DMRG method and have been 
compared with the corresponding spin-wave results. The spin of the ground state is given 
by sq = N(si — S2) for a 2N site system in all the three cases. There are two types of 
spin excitations in all the systems. The lowest excitation one is gapless in the infinite 
chain limit to a state with spin sq — 1. Another excitation to the state with spin sq + 1 
is gapped and the gap is larger when the difference in site spins \si — S2I is an integer. 
Interestingly, the low-energy spectrum remains qualitatively unchanged upon dimeriza- 
tion, and the dimerization is itself conditional in the infinite chain limit for all the three 
systems. 

We have also employed the DMRG technique to obtain low-temperature thermody- 
namic properties. For all these cases, the magnetic susceptibility shows very interesting 
magnetic field dependence. The yT vs. T plot shows a minimum at low-magnetic fields, 
and the minimum vanishes at high-magnetic fields. The temperature at which this min- 
imum occurs varies from system to system. The specific heat shows a maximum as a 
function of temperature at all fields. Moreover, the height of the maximum shows a 
dramatic increase at high-magnetic field. Interestingly, in each system, the temperature 
corresponding to the maximum in C and minimum in yT is the same at low-magnetic 
fields. 
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Figure Captions 



Figure 1 

Extrapolation of the ground state energy per site (e ), in units of J, as a function of 
inverse system size. 

Figure 2 

Expectation values of the z-components of the two spins vs. the unit cell index n. The 
upper and lower points are for the spin-sx and the spin-s 2 sites respectively. 

Figure 3 

Subtracted two-spin correlation functions as a function of distance between the two spins, 
(a) spin-sx spin-Sx correlations, (b) spin-s 2 spin-s 2 correlations, and (c) spin-sx spin-s 2 
correlations. In each figure, squares correspond to (|, 1)), circles to (|, |) and triangles 
to (1, \) systems. 

Figure 4 

Energy difference (units of J) between the ground state and the lowest energy state with 
spin s = sq — 1 as a function of inverse system size, sq is the total spin of the ground 
state. 

Figure 5 

Excitation gap (units of J) from the ground state (spin s = sq) to the state with spin 
s = sg + 1 as a function of the inverse system size. 

Figure 6 

Gain in magnetic energy AE(2N, 5) (units of J) associated with dimerization vs. the 
inverse system size for two different values of dimerization 5.5 — 0.025 (squares) and 
5 = 0.05 (triangles). 

Figure 7 

Magnetic energy gain (units of J) as a function of dimerization parameter 5 for different 
system sizes. In the figure 2N = 50 (squares), 2N = 100 (circles) and extrapolated values 
with N — > oo (triangles) are shown. 

Figure 8 

Log-log plot of extrapolated magnetic energy gain (units of J) for infinite system size and 
dimerization parameter 5. Interestingly, the slope is found to be 2.00 ± 0.01 for all the 
three systems. 

Figure 9 

Excitation gap (units of J) to the state with spin s = sq + 1 from the ground state 
(s = s G ) as a function of 5 for the dimerized alternating chain. The exponent is 1.0 it 0.01 
for all three systems. 

Figure 10 

Plot of magnetization per site as a function of temperature T for four different values of 
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the magnetic fields B. Squares are for B = O.lJ/gfis, circles for B = 0.5J/g/j,B, triangles 
for B = J/gfiB and diamonds for B = 2 J ' jg\iB- 
Figure 11 

Magnetization per site vs. the magnetic field strength B, in units of J/gfiB, for four 
different temperatures T. T = 0.3 J/ ks results are given by squares, T = 0.5 J/ ks by 
circles, T = 0.7 J jk B by triangles and T = J/k B by diamonds. 

Figure 12 

yT ( defined in the text) per site as a function of temperature T for various magnetic fields 
B. Zero field results are shown by squares, B = O.OU/gfiB by circles, B = O.lJ/gfiB by 
triangles and B = J/gfiB by diamonds. 

Figure 13 

Specific heat per site as a function of temperature T for four different values of magnetic 
fields B. Zero field data are shown by squares, B = 0.01J/gfx B by circles, B = 0.1J/gfj, B 
by triangles and B = J/g^B by diamonds. 
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